Dynamics of multifrequency oscillator communities 
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We consider a generalization of the Kuramoto model of coupled oscillators to the situation where 
communities of oscillators having essentially different natural frequencies interact. General equations 
describing possible resonances between the communities' frequencies are derived. The mostly simple 
situation of three resonantly interacting groups is analyzed in details. We find conditions for the 
mutual coupling to promote or suppress synchrony in individual populations, and present examples 
where interaction between communities leads to their synchrony, or to a partially asynchronous 
state, or to a chaotic dynamics of order parameters. 
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Networks of coupled oscillators describe synchroniza- 
tion in lasers [l| and Josephson junctions!!}, atomic 
recoil lasers Q, electrochemical oscillators [jj, applaud- 
ing persons in a large audience |5), pedestrians on foot- 
bridges .6] and many other systems. In the simplest 
setup, when all oscillators are of a same type and coupled 
via a mean field, the synchronization transition has been 
treated by Kuramoto [3, H| in analogy with the mean field 
theory of ferromagnetic phase transitions. Since that the 
Kuramoto model, where the coupled oscillators are rep- 
resented through the phase dynamics, has been used as 
a paradigmatic one for mutual synchronization of oscilla- 
tors [8l4lC|. In various generalizations this approach has 
been extended to the more complex and general situa- 
tions. One direction is introduction of complex coupling 
functions with possible nonlinear dependencies on 
the mean fields [l2-14|. Another very popular exten- 
sion of the Kuramoto model deals with heterogeneous 
oscillator populations, consisting of different communi- 
ties (groups) that differ in their contributions to the mean 
fields and in their response to these fields [l5l - [l8| . In par- 
ticular, nontrivial regimes appear if some interactions are 
"attractive" and other "repulsive" , or the oscillators can 
be characterized as "conformists" and "contrarians" . 

In most studies of the interacting oscillator communi- 
ties it is assumed that all oscillators have close frequencies 
around some basic one. For a small coupling this allows 
one to obtain, by virtue of averaging over the basic pe- 
riod, equations containing phase differences only, and to 
apply the Kuramoto method to them. In this letter we 
extend the theory to the case of multifrequency communi- 
ties, where natural frequencies of interacting groups differ 
significantly. In this case one cannot perform a common 
averaging, but has to check if resonances between differ- 
ent communities are present. In our previous study we 
have focused on the non- resonant case [2(J , in this letter 
we show that resonant interactions between communities 
lead to non-trivial effects of mutual synchronization and 
desynchronization of groups, and also to chaotic behavior 
of the mean fields. 



We start with a formulation of general phase equa- 
tions for resonantly interacting oscillator communities. 
Oscillators are described by their phases, and interact via 
mean fields, produced by communities. A field produced 
by a community with index m can be thus represented 
as a function of the generalized order parameters (2l| of 
this community 
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where () means averaging over the community members. 
A combination of these fields constitutes a force Q(t), 
acting on an oscillator from community 0, having fre- 
quency close to ll>q , which influences its phase dynamics 
as d©] 

<p = w + Aw + S(<p)Q(t) = w + Aw + ^ s n e mv Q(t) 

where S(<p) = s ne mv is the phase sensitivity function 
of the oscillator and Aw is a small individual deviation 
from loq- Representing this phase as ip = Lo$t + (p, where 
(p varies slowly on the time scale w^ 1 , we can average 
over the period 2-7r/wo to get 
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where q n u) = J Q " /u "° dt Q(t) exp[mwot] is the component 
of the forcing at the frequency — nujQ. In this component 
we have to consider the slowly varying ingredients of the 
forcing as the "frozen" ones. To separate slow and fast 
time scales in the forcing, we represent the order param- 
eters as 

_ 2,(rn) e iku m t 

where w m is a basic frequency of the community with 
index m, and Z are slow. Substituting this in Q and ex- 
panding in powers of order parameters, we can generally 
write 

<7n Wo = y^-Z| m) <5(fcw T , 
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We see that a direct interaction between communities m 
and is possible if their basic frequencies are in a rational 
relation kui m = uujq (if this relation is fulfilled only ap- 
proximately, one uses the freedom in the definition of the 
basic frequency of the community, and shifts u m slightly 
to have an exact resonance). Furthermore, the second 
term in ([3]) describes a "triplet" interaction if three com- 



munities have frequencies satisfying ku)„ 



nuj ; 



also higher-order interactions between four communities, 
described by cubic in Z terms, are possible, etc. 

In this letter we do not aim to consider all possible 
cases of resonant multifrequency interactions contained 
in (12l:j|) but focus on a simple example. We assume the 
simplest situation where the phase sensitivity function 
has only the basic harmonics n = ±1, and only the first- 
order mean fields with k = ±1 in ((TJ contribute to the 
coupling. In this case the basic resonant condition in- 
cludes three communities: wi + lu 2 = lu 3 ~ A, with small 
mismatch A. Taking also into account interactions in- 
side communities (which are described by the first term 
in ([3]) with u m = ujq), we end up with the phase model 
describing the resonant interaction of oscillators in three 
communities (cf. [22[ for a particular case of 1 : 2 reso- 
nance). To simplify notations, we denote the phases in 
these communities as <p k , ipk , Ok , and the corresponding 
order parameters as z\ 
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4> k = Ul + Awi, fc + 2Im[(e 1 z i + 7l z 2 *z 3 ) e - #fc ] 

ip k = 0J 2 + Aw 2 , fc + 2Im[(e 2 z 2 + liz\z 3 )e-^ k \ (4) 

9 k = uj 3 + Auj 3 . k + 2Im[(e 3 z 3 + j 3 z 1 z 2 )e~ iek } 

Here terms Awi_ 3i fc account for a distribution of frequen- 
cies of individual oscillators within communities, ej = 
Sie lai and 7$ = r^e 1 ^ are complex coupling constants. 
In the absence of mutual resonant coupling (r = 0), 
each community is described by the standard Kuramoto- 
Sakaguchi model [23j |. It is instructive to write micro- 
scopic phase equations for the phases: 

(j) k =ui + Auji ik + £i sm(<pj - <j) k + ai) + 
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ipk = u 2 + Aoj 2 ,k + £ 2 sin('0j - ipk + a 2)+ 
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in. I 

9 k = u 3 + Auj 3 , k + e 3 ^2 sin (°3 ~ °k + "3)+ 



T 3 2^ sin(0 m + - 9 k + ft) 



To obtain a closed system of equations for the or- 
der parameters zi.,2,3, we adopt the Ott-Antonsen ap- 
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proach 2J, |25[ , in which a particular form of the distri- 
bution of the phases is assumed, parameterized by the 
order parameter. If, furthermore, a Lorentzian distribu- 
tion of frequencies around the basic ones is considered 
(with widths #1,2.3)1 the Ott-Antonsen equations take es- 
pecially simple form: 

zi = zi(iut - Si) + (eizi + 71 z 2 z 3 - z\{t\z\ + 7l z 2 z 3 )) 

z 2 = z 2 (iu> 2 - 6 2 ) + (e 2 z 2 + 7 2 z^z 3 - z\{e\z% + 7 2 ziZg)) 

z 3 = z 3 (iu 3 - S 3 ) + (e 3 z 3 + 73 ziz 2 - z|(e 3 z 3 + 73^*^2)) 

(6) 

System ((6]) , describing the dynamics of order parameters 
of three resonantly interacting communities, is the main 
object of our analysis below. We focus on specific fea- 
tures resulting from the mutual coupling, where it acts 
"against" the internal coupling within the communities. 

Essential properties of the coupling between oscilla- 
tors such as their tendency to synchrony or to asyn- 
chrony, depend on the arguments of the complex coupling 
parameters oti and ft. For the coupling inside a com- 
munity, the argument a corresponds to the phase shift 
in the oscillator-to-oscillator coupling in the Kuramoto- 
Sakaguchi formulation [23j ]. for cos a > the coupling 
is attracting and synchronizing, while for cos a < it is 
repulsing and desynchronizing. A corresponding inter- 
pretation of arguments of mutual coupling ft is not so 
straightforward. To achieve it, we rewrite the complex 
system (|6]) in terms of the amplitudes and the phases of 
the complex order parameters Zj = pi exp[i$.;]: 

pi = S1P1 + (1 - Pi)(eipi cosai + T x p 2 p 3 cos(* + ft)) 
h = -5 2 p 2 + (1 - pl){e 2 p 2 cosa 2 + T 2 p 1 p 3 cos(^ - ft)) 
A3 = S 3 p 3 + (1 - p 3 )(£ 3 p 3 cosa 3 + T 3 pxp 2 cos(* - ft)) 
* = A' - (p^ 1 + p 3 )T 3PlP2 sin(* - ft)- 

(p^ 1 + p 2 )T 2PlP3 sin(* + ft)- 

(Pi 1 + Pi)T lP2 p 3 sin(#-ft) 

(7) 

Here A' = A + (1 + pDeasinas - (1 + p 2 )e 2 sina 2 - 
(1+ p\ )e% sinai is the effective frequency mismatch that 
includes also frequency shifts due to intra-communities 
interactions, * = $ 3 - $ 2 - $1 + 0.5(ft + ft) is the 
phase difference between the communities order param- 
eters, ft = 0.5(ft - ft) and ft = ft + 0.5(ft - ft). 
One can see from the equations for pi, that the effect of 
the inter-community coupling depends on signs of terms 
cos ('J ± ft), cos ('I' — ft). These depend on the dynam- 
ics of the phase difference so a general conclusion is 
hardly possible. Let us look on the simplest situation of 
exact resonance, where A' = 0, of equal coupling con- 
stants Ti = r 2 = r 3 , and of equal order parameters in 
each community p\ = p 2 = p 3 . Then the stable phase dif- 
ference is = arctan(sinft(cosft + 2cosft)~ 1 ). Sub- 
stituting this solution, we come to the diagram Fig. [TJ 
which shows the regions of positive and negative signs 
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FIG. 1. (Color online) Regions of synchronizing and desyn- 
chronizing effect from the triplet coupling, for different combi- 
nations of arguments of coupling constants /9i, /?2- The effect 
on a community is marked by (+) for enhancing synchrony, 
and by ( — ) for a desynchronizing action. In small overlaps of 
the ovals these markers "multiply" (so, e.g. in the overlap of 

ovals (H h) and (+ H — ) we have (H ), i.e. community 

1 synchronizes by coupling while communities 2,3 desynchro- 
nize). 



of factors cos^o ± ft), cos^o — ft), i.e. the regions 
where mutual resonance coupling between communi- 
ties promotes (for positive signs) synchrony or tends to 
desynchronize (for negative signs). In analytical form, 
these conditions are: Communities 1, 2, 3 synchronize for 
2 cos 2 ft + cos(ft + ft) > 0, 2 cos 2 ft + cos(ft - ft ) > 0, 
and 1 + 2 cos ft cos ft > 0, respectively, and desynchro- 
nize otherwise. We see that desynchronization in all three 
communities is not possible, while there are situations 
where two desynchronize, regimes where one desynchro- 
nizes, and states where mutual interaction improves syn- 
chrony in all communities. Below we give examples for 
synchronization and desynchronization effects in several 
setups. 

Mutual synchronization. In this part we assume that 
internally in each community the coupling is either re- 
pulsing or weakly (subcritically) attracting and without 
the mutual coupling the asynchronous states pi — are 
stable. To see that the mutual coupling can synchronize, 
we consider a simple symmetric case where the distri- 
bution widths and coupling constants are the same for 
all communities <5; = 5, Ei = e, Ti = T, the arguments 
of coupling constants vanish a-i = ft = 0, and there is 
no mismatch A = 0. Then it is easy to see from ([7]) 
that \ff — ► 0. Moreover, from the existence of a non- 
increasing on the trajectories of system (|7|) Lyapunov 

function L = -T Pl p 2 p3 - EiUd^M 1 ~ Pi) + £ Pi) [t 
follows that in this system only equilibria are possible. 
According to our assumption S > e, so the asynchronous 
state pi = p2 — P3 — is always stable, while for large 
enough T another synchronous state appears via a saddle- 
node bifurcation. In Fig. [2] we show the regions of pa- 
rameters with bistable synchrony-asynchrony states, and 
illustrate the appearance of the synchronous states as the 
mutual coupling is increased. 
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FIG. 2. (Color online) (a) Region on the plane of parameters 
where partially synchronous state appears, (b) Bifurcation 
diagrams showing dependence of the steady state order pa- 
rameters (equal for all communities) on the mutual coupling, 
for different couplings i — e/8 within the groups. 



Mutual desynchronization. As above, we assume here 
equal parameters for communities' heterogeneity 5 and 
the internal coupling e, so that the latter is real and ex- 
ceeds the critical value for synchronization e > S, also we 
set A = on — 0. To account for possibly desynchroniz- 
ing mutual interactions, we need to have non-zero argu- 
ments at least in some constants of mutual coupling. To 
simplify, we assume symmetry of the two low-frequency 
communities 1 and 2, setting T± — Y 2 — b, Tj, = c, 
and ft = 0. According to Fig. [TJ synchronization in 
community 3 may be destroyed if ft is close to tt. In 
order to clearly see the desynchronization effect of the 
mutual coupling, we assume S — > 0, in this case the syn- 
chronous communities are in fact identical clusters with 
order parameters p = 1, and a constant phase difference 
= — arctan(csinft(2& + ccosft)^ 1 ). This solution is 
stable for small ft, but loses stability through a trans- 
critical bifurcation at the critical value of this parame- 
ter: cosft, c = (1 - c 2 - V46c 2 - c 2 + l)(26c)- 1 . Beyond 
this transition, communities 1 and 2 remain in synchrony, 
while community 3 becomes partially synchronized, first 
with a constant order parameter, and at ft closer to 7r, 
with a periodically oscillating one. We show the bifurca- 
tion diagram for the partial desynchronization transition 
in Fig. H 

Chaotic order parameters. In the case when the mu- 
tual coupling between communities is much stronger than 
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FIG. 3. (Color online) Bifurcation diagrams showing depen- 
dence of the order parameter p3 on the coupling argument 
/?2 for Si = 1, c = 2, and different values of 6: 6 = 2 in (a) 
and b = 10 in (b). Solid and dashed lines are stable and un- 
stable equilibria; red line with markers denotes maxima and 
minima of periodic oscillations. TC, SN and AH denote a 
transcritical, a saddle-node, and an Andronov-Hopf bifurca- 
tions, respectively. 



the internal one, complex synchronization patterns in- 
cluding chaos are possible (for other examples of chaotic 
order parameters in coupled communities see [2614281 ] ) . 
We show in Fig.|4]an example of such a chaotic variation 
of the order parameters Pi,2,3 for the case where intra- 
communities couplings lead to synchrony in the groups 
e > S, but due to mutual interactions chaos appears in a 
certain range of arguments of mutual coupling /3i , /?2 ; for 
real values of the mutual coupling constants 7, i.e. for 
Pi = 0, we have not found complex behaviors. 

In conclusion, we have derived general equations de- 
scribing in the phase approximation the resonant inter- 
actions between communities of oscillators, which basic 
frequencies differ from each other, but are in a combi- 
national resonance. Such a situation can be observed in 
populations of neurons: it is known that in brain regu- 
lar macroscopic activity is observed across different fre- 
quency ranges 29]. Especially alpha, gamma, and theta 
bands may demonstrate quite regular oscillations, inter- 
action of which can be treated according to the presented 
framework. 

We focused in this letter on a detailed description 
of the most elementary three-community "triplet" reso- 
nance uj\ + u)2 ~ UI3, in terms of the evolution of commu- 
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FIG. 4. (a) Chaotic regime in system J7]) for /3i = 1.2. Time 
evolution of order parameters pi,2,3- (b) Results of calcula- 
tion of the largest Lyapunov exponent A for system J7J in 
dependence on parameter Nonzero values of A indicate 
chaotic regime. Parameters in both cases are the following: 
8i = 0.1, e, = 0.12, a, = 0, T t = 3, /3 2 = tt, A = 0.3. 



nities order parameters. This is accomplished by using 
the Ott-Antonsen ansatz allowing one to write a closed 
system for three complex order parameters. Remarkably, 
the inter-community interaction not only shifts relative 
phases of the communities mean holds, but influences in- 
ternal synchrony within communities. We have demon- 
strated how the inter-community interaction can induce 
or suppress internal synchronization. Furthermore, we 
have shown that resonant interaction of communities can 
lead to chaotic dynamics of the order parameters. 
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